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Aims. We investigate the circumstellar dust shell of the water fountain source IRAS 16342-3814. 

Methods. We performed two-dimensional radiative transfer modeling of the dust shell, taking into account previously 
observed spectral energy distributions (SEDs) and our new J-band imaging and H- and Ks-ha,nd imaging polarimetry 
obtained using the VLT/NACO instrument. 

Results. Previous observations expect an optically thick torus in the equatorial plane because of a striking bipolar 
appearance and a large viewing angle of 30 - 40°. However, models with such a torus as well as a bipolar lobe and an 
AGB shell cannot fit the SED and the images simultaneously. We find that an additional optically and geometrically 
thick disk located inside a massive torus solves this problem. The masses of the disk and the torus are estimated to be 
0.01 Mq at the a max = 100 /im dust and 1 Mq at a max = 10 /im dust, respectively. 

Conclusions. We discuss a possible formation scenario for the disk and torus based on a similar mechanism to the 
equatorial back flow. IRAS 16342-3814 is expected to undergo mass loss at a high rate. The radiation from the central 
star is shielded by the dust that was ejected in the subsequent mass loss event. As a result, the radiation pressure 
on dust particles cannot govern the motion of the particles anymore. The mass loss flow can be concentrated in the 
equatorial plane by help of an interaction, which might be the gravitational attraction by the companion, if it exists in 
IRAS 16342-3814. A fraction of the ejecta is captured in a circum-companion or circum-binary disk and the remains 
are escaping from the central star(s) and form the massive torus. 

Key words. Stars: AGB and post-AGB - circumstellar matter - radiative transfer - individual (IRAS 16342-3814) 



1. Introduction 

IRAS 16342-3814 (hereafter 116342), also known as OH 
344.1+5.8, is an oxygen-rich proto(pre)-planetar y nebula 
(PPN) exhibiting a striking pair of bipolar lobes. He et al.1 
pOOl detected a low 12 CO/ 13 CO intensity ratio of 1.7 
and concluded that 12 C is being converted into 13 C via 
the hot bottom burning process, suggesting that the initial 
mass of the central star is as high as ^5 Mq. Another 
characteristic of this object is high-velocity maser com- 
ponents. The velocities of OH masers at 1612, 1665, and 
1667 MHz are ^65 kms -1 , which are much faster than those 



in many OH/IR stars (jLikkel fc M orris! 1 1 9881). The H 2 Q 
masers are even faster ~ 130 - 180 kms -1 ( Zuckerman fc Lol 
ll987tlLikkel k Morrill 9881) . The apparent dynamical ages 
of the masers are very short at Si 100 yr (|Claussen et al.1 
20091 llmai et al.l l2ul2). The H2O masers are more like jets 
in appearance than the typical AGB winds. Because of 
these maser properties, o bjects like 116342 are called "water 
fountain source (WFS)" (|Likkel k MorridllQSallmai et al.1 
120071). So far, 14 objects are classified into this group 
(llmai et al.1 l2007t ISuarez et al.1 120091: IWalsh et al.l 120091: 
iGomez et al.ll2oillh Many PNs exhibit jets, which cannot 
be explained with onl y the generalised inter acting stellar 
wind (GISW) model dSahai k Traugerl 1 19981k The WFSs 
are thought to be progenitors that provide crucial clues to 
understanding the shaping mechanisms of these PNs. 

Because of the peculiar maser properties in WFS, many 
works have studied H2O and OH masers. In 116342, the 



best-studied object, the circumstellar dust shell (CDS) has 
been studied. The optical and near-infrared (NIR) im- 
ages show a distinct bipolar appearance wit h an exten- 
sion of ~5" lying at a position angle of 73° (jSahai et al.l 
Il999l ). In the optical, the eastern lobe is faint er by a fac- 
tor of ~0.1 with respect to the western one. ISahai et al.l 
(|l999l ) estimated the viewing angle to be ~40° from the 
1612 MHz OH maser distribution, i.e. the western po- 
lar axis is tilted by ~40° toward the observer. The Keck 
II i'-band im ages show a cork screw-like feature in the 
bipolar lobes (jSahai et al.1 120051 ) . Although a clear corre- 
lation with the distribution of the H2O maser jets is not 
found, the feature is probably a kind of precessing jet 
that carved out the inner part of the bipolar lobe. In the 
mid- infrared (MIR), a n elliptic flux structu re was detected 
(jMeixner et all 119991: iDiikstra et alll2003l ). However, the 
high-resolution VISIR/VLT data resolved a pair of bipolar 
lobes with a separation of the flux peaks of '/92 at 11.85 /im 
(jVerhoelst et al.l 120091: iLaeadec et~alll201lh . These results 
suggest an optically and geometrically thick torus in the 
equatorial plane of this object. 

In this present work, we have modeled the dust shell of 
116342 by means of two-dimensional radiative transfer cal- 
culations. The spectral energy distribution (SED) collected 
from various sources and our new NIR polarimetric data 
obtained using the VLT/NACO instrument were used to 
constrain the physical parameters of the CDS. In Sects. [2] 
and El our VLT/NACO observations and radiative transfer 
calculations are presented. In Sect.^J we discuss the struc- 
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ture and the formation of the inner most region of the dust 
shell. 



2. VLT/NACO imaging polarimetry 

2.1. Observations and data reduction 

We obtained standard (non-polarimetric) images in the J- 
band and polarimetric images in the H and K$ band using 
the CONICA camera with the NAOS adaptive optics (AO) 
system mounted on the Very Large Telescope (VLT). The 
S27 camera with a pixel scale of 27.15 mas pix -1 was oper- 
ated in the Double_RdRstRd read-out mode. For the wave- 
front sensing, the N20C80 dichroic filter was used. Because 
the target is extended in the optical and near-infrared and 
the central star feature is invisible, a single star 2MASS 
J16374030-3820087 (m K = 8.0 mag) at a 9"7 separation 
from the target was monitored for the AO reference. 

The J-band imaging was carried out on August 4, 2009. 
The natural seeing was 0'.'5 - 0'.'6 according to the ESO 
observatories ambient conditions data base. The exposure 
time was 30 sec. per frame and six images were taken at 
each four jitter positions, yielding a total integration time of 
12 min. The raw data were reduced by subtracting the dark 
frame, dividing by the flat frame, and combining the all 
science data. We used the photometric data of GSPC S875- 
C (to,/ = 11.085 mag) for the flux calibration, which was 
offered in the ESO standard calibration plan. The signal- 
to-noise ratios per beam is 40 - 240 in regions with the 
surface brightnesses of 2 x 10~ 15 Wm _2 /im _1 arcsec -2 or 
brighter. 

The H and K$ band imaging polarimetry was car- 
ried out on 2 May, 2009. The natural seeing varied be- 
tween l'.'O and l'/5 during the observations. The turnable 
half-waveplate (HWP) and the Wollaston prism with a slit 
width of 3" were used to measure linear polarization. The 
slit position is aligned along the long axis of the target at 
a position angle of 73°. Four image sets were taken at posi- 
tion angles of the HWP of 0°, 45°, 22.5°, and 67.5°. Seven 
and ten jitter positions were performed and exposure times 
were 60 sec. and 30 sec. per frame, in the H and K$ band, 
respectively. The total integration times were 28 min. in the 
H band and 20 min. in the K$ band. After the subtraction 
of the dark frames, the sky background levels and flat field- 
ing, the Stokes IQU parameter images were obtained by 
combining the ordinary and extra-ordinary ray frames. The 
data of a single star HD 147283 were used for polarization 
calibration (P H = .91 ± 0.05 % and P K = 0.48 ± 0.04 %; 
IWhittet et al.lli~992[ ). Our uncalibrated measurements are 
Ph = 3.6 % and Pr = 1-7 %. These discrepancies include 
the instrumental polarization and some other error sources. 
In observations that require a polarization accuracy better 
than 1 %, such as the lin e spectropolarimetr y, a careful cal- 
ibration is essential (see lWitzel et alj|2011|) . On the other 
hand, an ^2 % accuracy is acceptable in most cases of 
optical and infrared imaging polarimetry of bipolar reflec- 
tion nebulae. The difference of these tolerances also results 
from the different interests of the observations. In imag- 
ing polarimetry, the spatial information is important, but 
is affected by the photon noise and the temporal variations 
in the shape of the point-spread function (PSF) because of 
the use of small "aperture" sizes, e.g. beam sizes. It is oppo- 
site in the aperture polarimetry or line spectropolarimetry. 
We assumed that the discrepancies between our observa- 



Table 1. Parameters of the PSF model. 
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Fig. 1. Radial plot of the PSFs in the J (red), H (green), 
and Ks (blue) bands. The surface brightnesses are nor- 
malized. The observations are compared with the two one- 
dimensional PSF models of Moffat and Gaussian functions. 



tions and the previous aperture polarimetry (|Whittet et al.l 
Il992f ) are caused by the instrumental polarization and the 
offset Stokes Q/I and U/I values. The estimated signal- 
to-noise ratios per beam of the Stokes I images in regions 
with surface brightness of 2 x 10 -15 Wm -2 /im -1 arcsec -2 
or brighter are 60 - 120 in the H band and 20 - 120 in the 
K$ band. The errors of linear polarization are Pjj ~ 3 % 
and P K ~ 2 %. 

We evaluated the shapes and sizes of the beam. In the 
J band, a single star detected in the science frame was ex- 
amined. In the H and Ks band the standard star was used 
because no single stars are detected in the science frames. 
We found that most of flux from the reference stars is de- 
tected within l'/5 from the central star, and the star fea- 
ture is nearly circular. Therefore, the star feature was tan- 
gentially averaged and was examined as a one-dimensional 
PSF. Figured] shows the radial plot of the PSFs in the J, H 
and Ks bands. The solid curves are for the observations. 
The measured full width at half maximum (FWHM) are 
O'.'ll, 0'.'12, and 0'/07 in the J, H and K s bands, respec- 
tively. The real PSF sizes of the object would be somewhat 
larger than them because an offset AO ref erencing was ap- 
plied . We examined a Gauss ian model (e. g. lBendinelli et al.l 
I1990D and a Moffat model (lMoffatlll969h . The model PSF 



forms are given by Fbauss (?') = ^4exp 



-1/2 {r/*Y 



and 



Ft 



Moffat ' 



(r) = A/ l + {r/a)' 



, where a is the beam radius 

and f5 determines the overall shape of the PSF. The re- 
sulting parameters and the radial plot of the modeled PSF 
are listed in TableQ]and Fig-Hi respectively. The Gaussian 
function is only able to fit around the central star and can 
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Fig. 2. J, H, and Ks band images of IRAS 16342-3814 obtained using NACO on the VLT. The polarization vector lines 
are overlaid on the H and K$ band images. The field of view is 2'.'7 x 5". The images are rotated counterclockwise by 
107°. The color scale bar indicates the surface brightness in Wm -2 /im _1 arcsec -2 . 



provide a good estimate of the beam size. On the other 
hand, the Moffat function fits better throughout the entire 
PSF. Although the resulting model of this simple analysis 
would be inaccurate compared to more sophisticated im- 
age reconstruction methods such as the Richardson-Lucy 
deconvolution, it is useful for our purpose of comparing 
images produced by radiative transfer modeling with the 
observations. 



2.2. Results 

Figure [2] shows the intensity (Stokes /) images in the J, 
H, and K$ bands. The polarization vector lines are over- 
laid on the H- and Ks band images. The intensity images 
show a peanut-like pair of bipolar lobes with approximate 
extensions of 2"5 x 4" and the appearance is nearly con- 
stant in all bands. The western (upper) lobe is about three 
times brighter than the eastern (lower) one. T hese results 
are co nsistent with Keck II images presented bv lSahai et al.l 
(|2005f ) . Our images show uneven flux structures in the bipo- 
lar lobes but do not clearly detect the corkscrew-like struc- 
tures that are seen in the Keck II L p -band image. 

The degree of polarization reaches Ph ~ 25 % and 
P K ~ 20 % in the brighter lobe and P H ~ 38 % and 
Pk ~ 38 % in the fainter lobe. The polarization vectors 
are aligned along the equatorial plane. The physical rea- 
son is the PSF smoothing effect, namely, the highly scat- 
tered polarized light in the bipolar lobes spreads with an 
extended PSF halo component (see Fig.[T]). This cancels the 
polarization component perpendicular to the equatorial di- 
rection and the parallel component remains. Indeed, our 
model results convolved with the modeled PSF reproduce 
the vector alignment as seen in Fig.[S] 



One of the aims of our NIR imaging polarimetry was 
to detect the three-dimensional dusty structure in I16342's 
CDS. If the corkscrew-like pattern, as detected in Sahai's 
L p band image, is formed by a jet, the degree of polar- 
ization would gradually vary along the spiral because the 
scattering angle also varies along it. Another is the hollow 
structure in the bipolar lobe, which could be detected along 
the rim of the bipolar lobe as a polarization enhancement 
for the optically thin case or the opposite for the optically 
thick case. The hollowness is most likely the result of an an 
interaction between the fast post-AGB wind and th e slow 
AGB wind, as disc ussed in th e (G)ISW model (e.g . Kwokl 
Il98l iBalick et al.1 fl987t llckl[T985 iSoker fc Liviol 11989ft . 
Although an uneven distribution of polarization is detected 
in our data, we are not able to provide suitable interpre- 
tations that explain these factors. Therefore, we used the 
polarization data only to estimate the particle sizes in the 
bipolar lobes, but not to discuss the dust shell structure in 
the other sections. 



3. Dust shell model 

We have performed radiative transfer calculations of 
I16342 's CDS using our Monte Carlo code (jMurakawa et al.l 
l2008af l . The SED, and intensity and polari zation images are 
considered sim ultaneously, as done before (jMurakawa et al.l 
I2008bll2010ah . We first tried some test models to obtain ap- 
proximated solutions and to understand which parameters 
affect which results. Then, many parameter sets (several 
thousands in total) were examined by SED fit, which takes 
less computation time. From this result, a dozen good pa- 
rameter sets were chosen and their images were evaluated. 
We finally selected one that reproduced the shape of the 
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Table 2. Parameters of our radiative transfer modeling. 3.2. Model geometry 



parameters values 


comments 1 




central star 




T e a 


3000 K 


^ ooUU 


Li, 


6000 L Q 


2 


D 


2.0 kpc 


2 


R* 


2.0 x 10 13 cm 


calculated 




disk 




Rin 


20 


20 - 30 


-Rdisk 


200 AU 


100 - 200 


A/disk 


0.01 M Q 


adopted 3 


Q-max 


100.0 jum 


adopted 


TV 


410 


calculated 




torus 




-^torus 


1000 AU 


700 - 1000 


^^torus 


1.0 M Q 


0.7 - 1.0 4 


Q-max 


10.0 ^m 


adopted 


TV 


240 


calculated 




bipolar lobe and AGB shell 


-Rlobc 


5000 AU 


adopted 5 




3.0 


adopted 


7 


0.3 


adopted 


£in 


0.01 


adopted 




150 


100 - 200 


-Rout 


12 000 AU 


adopted 


M env 


1M Q 


1 - 1.5 




5.0 /jm 


2.0 - 5.0 



1 Ranges g ive the uncertainty of the corresponding model pa- 
rameters, 2 Sah ai et al.l (11999ft . 3 assuming an a max = 100.0 /im 
dust model, 4 assuming an a max = 10.0 pm dust model, 5 based 
on comparison of the intensity image. 



SED, bipolar appearance, and obscuration of the central 
star in the optical to NIR, and had nearly constant NIR 
polarizations in the bipolar lobe. 

Sect. 13. 11 13.21 and 13.31 describe the model assumptions 
and we compare the model results with the observations in 
Sect.EOl 



3.1. Stellar parameters 



15,000 K 
aperture 



The stellar temperature T c g was estimated to be 
from the optical photo metry with a 15" 
(jvan der Veen et aLI Il989h . whereas a strong 2.3 /im CO 
absorption was detected wi t h a 5 " aperture, indicating 
T off <, 3500 K. ISahai et al.l (|1999f l pointed out that this 
strong discrepancy is due to a contamination by the fluxes 
of two background stars (my =14. 9 and 13.9 mag) in the 15" 
aperture, which are brighter than 116342 (my=15.6 mag). 
Because the central star is invisible in the optical and in- 
frared, no precise investigation of the stellar properties have 
been reported. We assumed a T e g of 3000 K. The stellar lu- 
minosity and the distance are also not well constrained. 
Therefore, we adopted a luminosity of 6000 Lq, typical for 
AGB and post- AGB stars, and a di stance of 2 kpc, as con- 
sidered bv lSahai et~ai1 (|l999l . I2005D . The central star, as a 
radiation source, was assumed to have a blackbody spec- 
trum. We find that an effective temperature of 3000 K fits 
better than 15 000 K. 



Several geometry forms have been proposed for 
the CDSs around AGB and post-AGB stars (e.g. 

Meixner et Hi ll 9991: lUeta k Meixnerl 
120051). " 



Kah n fc West! Il98a 
2003| 



Qppcnhe imer et al.l 120051). Tho se proposed 

bv iKahn fc West! (I1985T T iMeixner et all (|1999D . and 



lUeta fc Meixnerl (|2003l ) have a latitudinal density gradient 
with a higher density at the equatorial region. Although 
these structures reproduce fan-shape morphologies for 
high equator-to-polar mass density ratios, e.g. ^10, the 
interaction of the stellar winds and an optically th ic k disk 
are not explicitly expressed. lOppenheimer et aD {2005) 
used one including an inner disk, a pair of bipolar lobes, 
and an outer spherical AGB shell. We employed a similar 
form for our previous work of a PPN M 1-92 where the 
bipolar appearance was reproduced well with this geometry 
(jMurakawa et al.ll2010"al) . 

We began with the latter form for 116342. However, 
we encountered a problem. The SED of this object has 
moderate NIR and strong MIR fluxes. At first glance, 
this (semi-) double-peaked flux structure can be explained 
with the standard cl a ssification of the post-AGB objects 
(|Hrivnak et al.l 119891: iKwokl I1993D . In this evolutionary 
phase, the intensive mass loss stops and the dust shell de- 
taches from the central star. This geometry with a large 
inner boundary causes less emission from hot dust at 3 to 
10 /*m and increases the scattered light in the optical-to- 
NIR. On the other hand, 116342 is expecte d to have a torus 
that is optically thick ev en in the MIR (jVerhoelst et al.l 
120091 : iLaeadec et al1l201ll ). With the aforementioned geom- 
etry form, we could not find a good solution to fit the SED 
and to reproduce a bipolar appearance in the NIR-to-MIR 
simultaneously. We attempted a geometry that includes an 
additional disk inside the inner disk, i.e. the inner region 
consists of an inner disk and outer torus. The mass density 
form of the entire dust shell is given by 



Pdisk (r, z) 
h(r) 

Ptorus 

(R) 

Plobc (R, 0) 
PAGE (R, d) 



Pdisk + Ptorus + Plobc + PAGB , 



Pa (r/ Rdisk) 2 25 exp 



2h 2 (r) J ' 



0.1i?disk (r/i?disk 
Pt (''/.Rdisk) -1 ' 8 , 

p c (R/R\ob c y 2 x 

p e (R/ R\obc) 2 , 



1.25 



^in ) 



(1) 

(2) 
(3) 

(4) 
(5) 



where the coordinates (R, 9) and (r, z) = (Rsin9, Rcos8) 
are the two-dimensional spherical and cylindrical coordi- 
nates. The disk and torus are placed in the regions where 
-Rin < R < -Rdisk and i? disk < R < -R t0 rus, respectively. The 
most plausible definition of a disk should be a region where 
the matter is gravitationally bounded by the central star 
and has a Keplerian rotating motion. Although dusty disk 
structures have been spatially resolved in severa l evolved 
stars (e.g. lMatsuura et al.ll2006t iLvkou et al.ll201ir) . the dy- 
namics are not well known or much less understood than 
accretion disks of young stellar objects. Therefore, we re- 
garded the disk simply as a geometrically thin, optically 
thick dust structure. In Sect. 13.41 we show that such an in- 
ner disk is required in 116342. 

The outer envelope consists of the bipolar lobe and 
spherical AGB shell. For the bipolar lobe, the regions 
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Fig. 3. Mass density map showing the inner 3000 AU region 
in a plane cutting through the symmetry axis. 

Rtoms < R < R\oh c (r {9) - 7 ) and R lohp {T (9) - 7) < 
R < -Riobc are the inner cavity and the rim of the lobe, 
respectively, w here T (9) = \9/tt — l/2f determines the 
lobe shape fsee iQppenheimer et al J 12003 ). The AGB shell 
is placed outside the bipolar lobe with an outer boundary of 
R ou t. The mass density coefficients are determined by the 
masses of the corresponding components. The dust mass is 
converted by assuming a gas-to-dust mass ratio of 200. 

Figure[3] shows the mass density map in the inner 
3000x3000 AU region. Our modeled geometry is compli- 
cated and it is in essence impossible to determine all pa- 
rameters accurately from our observational data. Therefore, 
some parameters are set to be fixed values. Important pa- 
rameters such as the inner disk radius i?disk> the disk mass 
Mjisk, the torus radius i?torus, and the torus mass M torus 
are estimated from our radiative transfer calculations. 

3.3. Dust particles 

iDiikstra et al.l (|2003l) discussed the mineralogy of 116342 
based on their ISO spectra obtained using the Short 
Wavelength Spectrometer (SWS). Their data show a vari- 
ety of absorption and emission features attributed to amor- 
phous and crystalline silicate, frostetite, crystalline water 
ice, etc. Investigating the dust chemistry is one of the in- 
teresting and important tasks in dust formation of evolved 
stars. However, a thorough treatment of so many differ- 
ent chemical compositions and identification of the spa- 
tial distributions of individual species are very compli- 
cated tasks and are, therefore, beyond the scope of this 
paper. In our modeling, we adopted simple dust models. 
We used the optical co nstants for oxygen-deficient silicate 
(|Ossenkopf et all [1991 . 

We assumed different particle sizes for the outer enve- 
lope, the torus, and the inner disk, but the same sizes for 
the bipolar lobes and AGB shell. The particle size in the 
outer envelope was estimated with the polarization data. 
Since the polarizations are high and do not change much in 



the H- and -fCs-bands, the fraction of small particles is ex- 
pected to be higher than that of the interstellar population. 
In our separated paper (|Murakawa et al.ll2012f) . we found 
a similar result in PPN 118276 and that a size distribution 
function with a steep power index provides a good estimate: 
0.05 fim. < a and n (a) oc a~ 5 ' 5 exp (— a/a c ), where a c is the 
cut-off size (|Kim et al.lll994l) . In 116342, a c = 5.0 /an fits 
the H and K$ band polarizations. 

Although large grains are expected in the disk and 
torus, we do not have (sub-)millimeter flux data or their 
gas masses and therefore we cannot constrain the sizes. 
Therefore, we provisionally adopted a size distribution 
of 0.005 ami < a < q ma x Mm and n (a) oc a~ 3 - 5 
(jMathis. Rumple fc Nordsiecldll977t ). a max for the disk and 
torus are set to be 100 /xm and 10 /mi, respectively. We will 
not discuss dust growth in this object with our results. 

3.4. Results 

Figure^] shows the SEDs. The black lines are the result of 
the selected model (model 1). The viewing angle was cho- 
sen to be 30° measured from the equatorial plane. Although 
this result is slightly different compared to that (40°) es- 
timated bv lSahai et al.l (|l999t) . this discrepancy is not im- 
portant in the following discussion. The model result fits 
the observation well from the optical to FIR, except for the 
10 /xm silicate absorption feature. 

We justify the presence of a geometrically thin, inner 
disk. First, we compared with a model without an inner 
disk (model 2), in which Mdisk is set to be and the other 
parameters arc identical. The difference is obvious in the 
wavelength range between 3 /xm and 10 /im. Although the 
disk mass in model 1 is only 0.01 Mq, the MIR flux is 
governed by the thermal emission from the hot dust in the 
inner disk. The other consideration is a model without an 
inner disk where we tried to fit the observations as best 
we can (model 3). This is our first attempt in the modeling 
(see Sect. l3.2[) . We find that models with small inner radii (a 
few tens AU) show too strong thermal emission at 3 /xm to 
10 /xm and that a model with a large radius of R[ n = 500 AU 
and Mtorus = 2.0 Mq fits best. However, as Fig. 2] shows, 
the flux between 3 /im to 6 /im is too low and model 1 
fits better. The interpretation of I16342's inner region is 
summarized as follows. An optically and geometrically thick 
torus should exist, as previous observations have shown. 
The torus is expected to have a large inner radius (a few 
hundred AU) and a high mass (1 Mq). An additional inner 
geometrically thin disk is required, which is responsible for 
the thermal emission in the 3 /im to 10 /xm. 

Figure[5] shows the monochromatic model images at 
1.25 /jm (J-band), 1.65 /im (iJ-band), and 2.05 /xm (K$- 
band) from left to right, which are simulated by scattering 
and absorption only. In these images, the Moffat PSF model 
images (see Sect. l2.H are convolved. We confirmed that 
the Gaussian PSF convolution produces a centro-symmetric 
polarization pattern in the bipolar lobe, but the vectors 
are aligned with the Moffat PSF model. This is quite un- 
derstandable because the surface brightness of a Gaussian 
function decays much faster than a Moffat function. The 
vector alignment in the entire nebula is reproduced by the 
PSF smoothing effect. Compared with the observations, the 
striking bipolar appearance and the presence of a dark lane 
in the equatorial plane are well reproduced. Although we do 
not present the results, we also confirmed that our model 



5 



Murakawa et al.: dust shell model of 116342 




1 ° / V 

wavelength \firr\) 








Fig. 4. Comparison of SEDs. Black 
lines are the result of the model with 
an inner disk. The solid, dotted, and 
dash-dot lines denote the total flux, the 
scattered light, and the thermal emis- 
sion, respectively. Red and blue lines 
denote the results of the same geome- 
try form but the disk mass of and a 
model with only a torus with an inner 
radius of 500 AU. Colored dots are col- 
lected from previous observations: US 
Naval Observatory (USNO), the HST 
Guide Star Catalog (GSC)-II, 2MASS 
all-sky catalog of point source, Catalog 
of Infrared Obser vations, Edition 5 
(|Gezari et all 1 19991 ). ISO photometry, 
SWS and LWS, MSX infrared point 
source catalog, IRAS point source cat- 
alog, a nd Akari/IRC mid- IR all-sky 
Survey (|Ishihara et al.ll2010f ) 



n 




offset (arcsec) 



offset (arcsec) 



offset (arcsec) 



Fig. 5. Results of our two-dimensional radiative transfer calculations. The wavelengths are the 1.25 /jm (J-band), 
1.65 /xm (ii-band), and 2.05 /im (iCg-band) from the left panels. The color scale bar indicate the surface brightness 
in Wm _2 ^m _1 arcsec -2 . The Moffat PSF model images are convolved (see Sect. 12. 1[) . 



image in the 11. 85 ^im shows a b i polar appe arance, as seen 
in the results of IVerhoelst et al.l (|2009h and iLaeadec et all 

poll . 



4. Discussion 

The mechanisms of PN shaping have been discussed in a 
number of publications before. To form a striking bipolar- 
ity with an equatorial narrow waist, the mass loss mate- 



rial should be highly concentrated in the equatorial plane. 
Morris (1987) modeled mass loss in a binary system from 
the primary red giant. Depending on the physical condi- 
tions, e.g. the mass and the separation of the secondary, 
a part of the mass loss material from the primary can 
be transferred onto the secondary, which forms a circum- 
companion disk. A part of accreting matter can be blown 
along the polar direction as a jet. Many bipolar PNs and 
PPNs, presumably 116342 too, are thought to form by the 
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binary interaction based on this idea. Both theories and ob- 
servations have developed the binary interaction hypothe- 
sis. 

Silicate carbon stars are objects in which oxygen-rich 
species exist around carbon stars. According to the cur- 
rent interpretation of the dual chemistry, the silicate dust 
that was ejected when the central star was an M giant is 
stored in a disk around the hypothetical binary companion 
and then the central star is t ransformed into a carbon star 
through the third dredge-up (|Llovd-Evanslll990ft . V778 Cyg 
is an object of this cl ass, which shows the 10 /im silicate 
emission feature, see lYamamura et all (|2000f ). These au- 
thors found that the dust has sub-micron sizes and temper- 
atures of 300 - 600 K and is located at about 24 AU from 
the primary. The present-day mass-loss rate is estimated 
to be about 10~ 8 MQyr -1 . They interpret that the mass 
loss material can be trapped in the circum-companion disk, 
but leaves soon after within a few orbital periods because 
the disk is optically thin and the dust is blown away by 
the strong radiation pressure from the primary. With these 
low mass-loss rates, the disks are unlikely to be responsi- 
ble for a bipolar appearance. It is of interest to mention 
the results of the smoothed particle hydrodynamic simu- 
lat ions for high mass-loss rate cases o f 10 -5 Mgyr" 1 
bv lMastrodemos k, Morris! (|l998l . ll999l hereafter M98 and 
M99). Their simulations assumed a mass-losing AGB pri- 
mary and a main-sequence companion (Ms = 0.25 - 2 Mq) 
with orbital separations of 3.6 - 50 AU. In M 1, 2, 10, 12, 
and 15 out of their 18 models, the density ratios of the 
equator to the pole Peq/Ppoi exceeds ~10 and the nebulae 
are classified as bipolar. M 1 and 2 are extreme cases, where 
the separations of the companion is small (6.3 AU in their 
models) and ratios of the accretion rate to the mass-loss 
rate are as high as ~0.1, yielding M acc *b 10~ 6 MQyr . 
Peq/Ppoi exceeds 100. These systems are likely to be pro- 
genitors of bipolar PNs with narrow waists. The thick 
winds sufficiently shield the strong radiation from the pri- 
mary and the mass loss matter is stabilized there for a 
long time (|Yamamura et al. 20001: ISo ker 2000 ), as expected 
in th e Red Rectangle (| Jura et al.l 1995L Ivan Winckel et alJ 
Il998j) . AFGL 2688 (e.g! iBieeing fc Nguven-Q-Rieul 11996ft . 
and probably 116342. 

We now consider a simple model of AGB mass loss to see 
the competition of the radiation and the shielding based on 
a similar co ncept to "equatorial back-flow" introduced by 
iSokerl (|2000f) . We assumed that a dust particle is directly ir- 
radiated by the radiation from the central star. The ratio of 
the radiation pressure to the gravitational attraction /3d us t 
is given by /3dust = L+Crp/AirGM+c, where C rp is the cross 
section of the radiation pressure per unit mass, G the grav- 
ity constant and c the velocity of light. Applying the stellar 
parameters of = 6000 Lq and M± = 1 Mq , a dust model 
that is derived from the outer envelope in our modeling, 
C rp =6.73x 10 3 cm 2 g _1 , where the stellar flux is assumed 
to dominate at A = 1 fim, /?dust is found to be ^3xl0 3 . 
This high value is quite reasonable in AGB winds, where 
the dust is blown away by the radiation pressure. Next, the 
shielding by the mass loss ejecta is considered. We assumed 
a spherically symmetric mass loss at a constant M and con- 
stant radial expansion velocity V r . The mass density has a 
power law distribution p (r) cx r -2 , where r is the distance 
from the central star. The optical depth between the inner 



boundary (r = R lu ) and a position at r = R is given by 



a*tM ( i 



47rV r 7 V^i 



(6) 



where C cx t is the extinction cross section per unit mass of 
the dust and 7 the gas-to-dust mass ratio of 200. If the 
mass-loss rate reaches 10 -5 Mgyr -1 at an expansion ve- 
locity of 15 kms -1 , and R 3> R ln = 5 AU is assumed, we 
derive an optical depth 



! 20.3 



15 kms 1 

V P 



M 



10- 5 MQyr- 1 



'5 AU 



(7) 



where C oxt = 9.10xl0 3 cm 2 g~ 1 . The dust shell is opti- 
cally thick. When dust particles move away, the radiation 
is shielded by the dust that is formed by the subsequent 
mass loss event. The optical depth attains ~8.0(= ln/3d U st) 
at R ~ 8.3 AU and the attenuated /3dust values becomes 
~1. The radiation pressure force cannot govern the motion 
of the particles. If the system has a close binary, the grav- 
itational attraction of the secondary can affect the motion 
of the mass loss ejecta and can concentrate the flow toward 
the orbital plane. This leads to a flatter geometry, i.e. disk. 

We have a rough scenario of the inner part of I16342's 
CDS as follows. The inner disk of 116342 in our model, 
which is introduced to be responsible for the MIR flux, is a 
circum-companion or circum-binary disk that is formed by 
this mechanism. Because the optical depth of the disk is so 
high at ^400 in the U-band that the effect of the radiation 
is negligible, the disk can be gravitationally bounded by the 
central star(s). On the other hand, the torus is formed in a 
somewhat different way. If we assumed that 116342 has a bi- 
nary companion, an angular momentum that is supplied by 
a close companion (say Ms ~ 1 Mq) is insufficient for the 
entire torus to have a keplerian rotating motion. Therefore, 
the torus is probably a structure that is accumulated from 
dust that is escaping from the central star(s). 



5. Conclusion 

We modeled the CDS of the WFS IRAS 16342-3814 with 
two-dimensional radiative transfer calculations. The model 
parameters were constrained by fitting the SED collected 
from various sources and the polarization images obtained 
using the VLT/NACO instrument. 116342 shows a striking 
bipolar appearance from the optical to the MIR wavelength 
ranges. Taking into account a large viewing angle of 30° - 
40°, an optically thick dust torus structure is expected in 
the equatorial plane. We first tried to fit a model consist- 
ing of a geometrically thick torus and an outer envelope 
with bipolar lobes and a spherical AGB shell. However, 
these models cannot reproduce the SED and the bipolar 
appearance simultaneously. We found that an additional 
component of an optically thick, geometrically thin disk is 
required inside the torus. In our modeling, different dust 
sizes were assumed for the disk, the torus, and the outer 
envelope. The particle sizes in the outer envelope were esti- 
mated from the polarization data. The polarization reaches 
25 - 38 % in the H band and 20 - 38 % in the K s band 
in our observations. This weak wavelength dependence is 
explained with a steeper power index of —5.5 of the size 
distribution function than that of —3.5 for the ISM. Similar 
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results are found in some other PPNs. Because of the sub- 
millimcter and millimeter flux information and the esti- 
mation of the disk are missing, the particle sizes of these 
component cannot be constrained well. We adopted dust 
models of a max = 100 /zm for the disk and a max = 10 fim 
for the torus. With this assumption, the masses of these 
components were determined to be 0.01 Mq and 1 Mq, 
respectively, but have large uncertainties. 

We discussed a possible formation scenario for the disk 
and torus. The important factors are the mass-loss rate and 
the binary interaction. For low- to intermediate-rates, the 
system is likely to have an optically thin circum-companion 
disk. To form an optically thick disk, the mass-loss rate 
would need to be so high that the thick wind blocks the 
radiation from the primary. This allows a part of the mass 
loss ejecta to be condensed in the equatorial region and to 
stay there for a long time. The remains, which are leaving 
from the central star(s), form the torus. 
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